Chaos, entanglement and decoherence in the quantum kicked top 
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We analyze the interplay of chaos, entanglement and decoherence in a system of qubits whose 
collective behaviour is that of a quantum kicked top. The dynamical entanglement between a 
single qubit and the rest can be calculated from the mean of the collective spin operators. This 
allows the possibility of efficiently measuring entanglement dynamics in an experimental setting. We 
consider a deeply quantum regime and show that signatures of chaos are present in the dynamical 
entanglement for parameters accessible in an experiment that we propose using cold atoms. The 
evolution of the entanglement depends on the support of the initial state on regular versus chaotic 
Floquet eigenstates, whose phase-space distributions are concentrated on the corresponding regular 
or chaotic eigenstructures. We include the effect of decoherence via a realistic model and show that 
the signatures of chaos in the entanglement dynamics persist in the presence of decoherence. In 
addition, the classical chaos affects the decoherence rate itself. 

PACS numbers: 05.45.Mt, 03.67.Bg, 42.50.Ct, 03.65.Yz, 32.80.Qk 



I. INTRODUCTION 

In classical mechanics, the chaotic behavior predicted 
for non-integrable systems can be qualitatively different 
from the regular dynamics of integrable systems. The 
concept of regular versus chaotic dynamics at the quan- 
tum level has been more difficult to define as there is 
no clear measure of chaos in the quantum regime. Un- 
derstanding the correspondence between quantum and 
classical evolutions in chaotic systems is a central prob- 
lem in quantum mechanics and a major focus of the 
field of quantum chaos. From a practical perspective, 
recent work [T] has shown that classical chaos can affect 
the implementation of quantum computing algorithms, 
and has fueled interest in identifying the effects of chaos 
on quantum information theoretic properties such as en- 
tanglement and fidelity. Entanglement is thought to be 
a fundamental resource for many quantum information 
processing applications, and the effect of chaos on the 
dynamical generation of entanglement has been a topic 
of several studies El El SI El IE1 IZ1 El (HI EU1 . The presence 
of chaos can also increase the rate of entanglement gener- 
ation between a system and its environment |12j , leading 
to increased decoherence and possibly stricter limitations 
on coherent quantum information processing. 

In this paper, we explore the effect of chaos on en- 
tanglement and decoherence in a quantum kicked top. 
We show that entanglement can be efficiently measured 
in this system, and identify signatures of chaos in the 
entanglement dynamics in a deeply quantum regime. 
The quantum kicked top [T31 [2] has become a stan- 
dard paradigm for theoretical studies of quantum chaos 
but has not yet been studied in experiments. Here, we 
propose and analyze a possible experimental realization 
based on Cesium atomic spins interacting with laser light 



and a pulsed magnetic field. With a ground hyperfine 
spin of F — 4 this system lies far from the semiclassical 
regime that is usually considered, since the size of h rel- 
ative to the total phase space is roughly 0.1. We study 
entanglement in the kicked top for parameters accessi- 
ble with this system. Our analysis shows that signatures 
of chaos can be observed even if the system undergoes 
decoherence and that chaos affects the decoherence rate 
itself. 

A quantum kicked top with total angular momentum 
j can be decomposed into N = 2j spin-1/2 subsystems 
(qubits) [6]. We consider here the entanglement between 
a single qubit and the remaining k — N—l qubits (hence- 
forth called 1 : k entanglement). The 1 : k entangle- 
ment is of relevance to our proposed experiment with cold 
atoms because, as we show, it corresponds to entangle- 
ment between electron and nuclear spin in a single atom. 
Measurement of entanglement is usually experimentally 
challenging as it would require performing complete state 
tomography. For the case of the kicked top, we derive a 
simple expression for a 1 : k entanglement measure in 
terms of the expectation values of the total spin of the 
N qubits. Hence, the 1 : k entanglement can be effi- 
ciently monitored by simply measuring the evolution of 
the mean total spin vector, thereby avoiding the need 
to perform complete state tomography. Although this 
result is generally applicable to any kicked top experi- 
mental realization, it is particularly useful for the cold 
atom experiments proposed here, in which the evolution 
of the mean spin vector can be efficiently measured in 
real time [15| 116]. 

We have performed numerical simulations of the kicked 
top dynamics to study the evolution of entanglement. 
Our results show that the long term 1 : k entanglement 
dynamics exhibit signatures of classical chaos even in the 
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extreme quantum regime considered here. Initial states 
localized in regular regions of the classical phase space 
exhibit quasiperiodic collapses and revivals, whereas ini- 
tial localized states centered in chaotic regions quickly 
spread out and do not exhibit quasiperiodic entangle- 
ment dynamics. These signatures are similar to those 
previously identified in the 2 : k entanglement [B] be- 
tween 2 qubits and the rest in a semi-classical regime. 
We explain the entanglement behavior by extending the 
analysis in [7] to the quantum regime. Quasiperiodic- 
ity versus irrregular dynamics depends on the support of 
the initial state on regular versus chaotic Floquet eigen- 
states as shown in [f\. The differences in eigenstate de- 
composition for initial states centered in regular versus 
chaotic regions are limited by the small size of the spin 
considered here, and has no significant effect on the ini- 
tial rate of entanglement generation as is the case for the 
semiclasscial regime 7] ■ Nevertheless we find surprisingly 
clear signatures of regular versus chaotic dynamics in the 
long term evolutions. Regular and chaotic structures of 
the mixed classical phase space are also clearly evident 
in the Husimi distributions of the Floquet eigenstates in 
the quantum regime. 

In addition to unitary evolution, we also consider the 
more realistic case of non-unitary dynamics when the 
kicked top system is coupled to its environment. Our 
simulations of the entanglement dynamics include the ef- 
fect of decoherence through a master equation descrip- 
tion that accurately simulates photon scattering in our 
cold atom system. Analysis of the negativity, a measure 
of entanglement for mixed states, shows that, although 
decoherence due to entanglement with the environment 
acts to reduce the overall 1 : k entanglement, striking dif- 
ferences in the entanglement dynamics in regular versus 
chaotic regimes can persist for times longer than the de- 
coherence time. Furthermore, the rate of decoherence it- 
self is slower in a regular regime than in a chaotic regime. 
Signatures of chaos in the decoherence rate were previ- 
ously analyzed in 112] using a generic model of deco- 
herence in a semiclassical regime. Here we use a realistic 
and accurate model to show that chaos affects the deco- 
herence rate in a physical system in a deeply quantum 
regime. Our proposed cold atom implementation would 
be the first to allow studies of decoherence and entangle- 
ment dynamics in a chaotic system. 

The paper is organized as follows: In Section II we de- 
scribe the basic features of the standard quantum kicked 
top and the corresponding classical system, which can 
exhibit chaos for certain dynamical parameters. Section 
III discusses the efficient measurement of 1 : k entangle- 
ment dynamics and the possibility of experimental stud- 
ies of entanglement with cold atoms. In Section IV we 
present an analysis of signatures of chaos in the entangle- 
ment dynamics of the kicked top for a regime accessible 
in the cold atom system. We build on our previous anal- 
ysis of the kicked top j7j and identify regular and chaotic 
eigenstates whose distributions are concentrated on the 
classical phase space structures in Section V. The rela- 



tionship between decoherence, entanglement and chaos is 
analyzed in Section VI. We present a summary and our 
conclusions in Section VII. 



II. THE QUANTUM KICKED TOP AND ITS 
CLASSICAL LIMIT 



The Hamiltonian for a quantum kicked top is given 

by nam 

oo 

H=^-J 2 x +pJ y J2 S{t-nr). (1) 

n— — oo 

Here, the operators J Xl J y and J z are angular momentum 
operators obeying the commutation relation [Ji,Jj] — 
iheijkJk- The Hamiltonian describes a series of kicks 
given by the linear J y term interspersed with torsions 
due to the nonlinear J, 2 term. The time between kicks 
is t, the angle of turn per kick is given by p and the 
strength of the twist is determined by k. The magnitude 
J 2 = j(j + l)h 2 is a constant of the motion 

The classical map from kick to kick can be obtained 
from the Heisenberg evolution equations for the expecta- 
tion values of the angular momentum operators, 

(Ji)n+1 = (U*JiU)n. (2) 

U is the Floquet operator describing unitary evolution 
from kick to kick, 

U = cxp(-iKj 2 /2j)exp(-ipJ y ), (3) 

where the energy is henceforth rescaled such that r = 1 . 
The equations describing the evolution of the mean J x , 
J y and J z involve functions of second moments of the 
angular momentum operators, whose evolution in turn 
depends on third moments in an infinite hierarchy of 
equations. In order to obtain the classical mapping, we 
factorize all second and higher moments into products of 
the mean values (1st moments) of the angular momen- 
tum operators. This corresponds to localization of the 
state to a point in the classical limit and breaks the infi- 
nite hierarchy of equations at the level of first moments. 
Then, by defining the normalized variables X — (J x )/j, 
Y = (J y )/j and Z = (J z )/j, we can write down the 
classical mapping 



X n +i 


= X n cos(p) + Z n sin(p) 


Y n+1 


= Y n cos(nX n ) - Z n sin(nX n ) 


Z n +i 


= Z n cos(KX n ) + Y n sin(KX n ) 


Z n 


= Z n cos(p) - X n sin(p) 


x n 


= X n cos(p) + Z n sm(p). 



X n and Z n are the angular momentum variables subse- 
quent to the kick but before the action of the twist. 
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FIG. 1: Stroboscopic phase space maps at different values of k, 
for kick strength p — n/2. The classical spherical coordinates 
(0, 4>) are plotted after each kick for 144 initial conditions, 
each evolved for 150 kicks, (a) For smaller k, the phase space 
is dominated by regular orbits. (b)At larger k, regular orbits 
are embedded in a sea of chaotic trajectories. 



The parameter n is the chaoticity parameter in the 
classical kicked top. For p = n/2, Fig. 1 shows how n af- 
fects the kick-to-kick stroboscopic dynamics of the classi- 
cal variables 8 — cos _1 Z, <f> — tan (Y/X). When k = 1, 
the stroboscopic phase space is dominated by islands of 
regular (periodic) motion [Fig. 1(a)]. As k is increased, 
the phase space becomes mixed with regular islands em- 
bedded in a sea of chaotic trajectories [Fig 1(b)]. For 
even larger values of k, the system eventually becomes 
globally chaotic. 



III. MEASUREMENT OF ENTANGLEMENT IN 
THE KICKED TOP 



where a ai are the Pauli operators for the ith qubit. For 
states in the symmetric subspace, the kicked top operator 
acts collectively on all N qubits preserving the symmetry 
of the N qubit state. This symmetry allows us to write 
the spin expectation values for any single qubit as 
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<Ja) 
2j 



(6) 



For an overall pure state, the entropy of the reduced 
state p, of a single qubit is a measure of 1 : k entangle- 
ment between a single qubit and the remaining k = N — 1 
qubits. For convenience we study the linear entropy, 
5=1 — Tr[/5 2 ], which ranges from for separable states 
to 1/2 for maximally entangled states. The single qubit 
reduced state p can in general be expressed as 



P 



(7) 



where s is the mean spin vector of the qubit. Using 
Eq (6), p can be expressed in terms of the collective an- 
gular momentum operators as 



1 



1 



P = 



cr. 



2 2j 

The linear entropy is then easily computed to be 



S = 



1 



1 



i-^((j x y + (j y y + (j z y) 



(8) 



(9) 



Measurement of the ensemble averaged angular mo- 
menta, (J Q ), thus provides us with enough information 
to calculate S, which is a measure of the 1 : k entan- 
glement, without having to perform complete state to- 
mography to reconstruct the entire multiqubit density 
operator. We also note that the function S is identical to 
the generalized entanglement with respect to the angular 
momentum observables - an entanglement measure that 
is independent of subsystem division |17j . Furthermore, 
the linear entropy function is directly related to the 'ex- 
tent' or average spread of the state on the sphere. For 
a state that is highly localized on the sphere, the sum 
of the mean values of the angular momentum operators 
with be close to j so the entropy S will be close to 0. On 
the other hand, for a highly delocalized state, the mean 
values of the angular momentum operators in Eq. (9) ap- 
proaches and the value of S approaches its maximum 
of 1/2. 



A. l:k entanglement measure 



B. Experiments with cold atoms 



We now focus on the dynamical evolution of entangle- 
ment in the quantum kicked top. The total spin j can be 
considered as a system of N = 2j qubits with 



N 

j _ \ " a otj 
i=l 



(5) 



An attractive physical system in which to realize 
a quantum kicked top is the total hyperfine (elec- 
tron+nuclear) F of the electronic ground state of an indi- 
vidual atom. In particular, we propose using samples of 
laser cooled alkali atoms to perform experiments on en- 
sembles of identical kicked tops. The atomic spins can be 



4 



initialized by optical pumping, and manipulated in a con- 
trolled manner by Larmor precession in applied magnetic 
fields, and electric dipole interaction with an applied laser 
field. Furthermore, one can perform polarization spec- 
troscopy on the driving field to probe the spins in real 
time without perturbing the dynamics [15]. These tools 
have already been used to demonstrate non-linear spin 
dynamics quantum state reconstruction |18j and 

quantum state control |19j with Cesium atomic spins. 
In the following, we briefly describe how to realize the 
kicked top with this system. 

In the low saturation, large detuning limit, the electric 
dipole interaction between a single atom and a monochro- 
matic laser field is described by a light shift 

U = -~E* ■ a ■ E, (10) 

where E = Rc(Ee~ lult ) is the electric field, and where 
the atomic polarizability a is a rank-two tensor operator 
that can be decomposed into irreducible components of 
rank 0, 1 and 2 [201 El] ■ Here, we consider alkali atoms 
restricted to a hyperfine ground state of given F. The 
light shift U is then an operator acting in a 2F + 1- 
dimensional manifold, separable into three contributions 
from the irreducible components of a. 

The rank-0 contribution is a scalar interaction which 
does not couple to the spin degrees of freedom and there- 
fore can be ignored. The rank-1 contribution is an effec- 
tive Zeeman interaction of the form B e g • F, where _B e ff 
is proportional to the cllipticity of the laser field polar- 
ization. For a linearly polarized driving field, this term 
disappears, leaving only the rank-2 contribution. Choos- 
ing linear polarization in the x-direction the overall light 
shift is then quadratic in a component of the hyperfine 
spin as required for the kicked top, 

U = f3h ls Fl (11) 

Here, 7 S = sr/2 is the single-atom photon scattering 
rate, with the saturation parameter s depending on the 
laser intensity /, detuning A , transition linewidth T 
and saturation intensity Joass = (r/2A) 2 (///o)- The 
parameter (3 is a measure of the relative timescales for 
unitary evolution and decoherence, and depends on the 
atomic species and the frequency of the driving field. 
It takes on a maximum value of 8.2 for Cesium atoms 
driven at a frequency in-between the two hyperfine com- 
ponents of the Dl line at 894 nm [TS]. Comparing Eqs (1) 
and(ll), we see that the strength of the chaoticity param- 
eter is related to the system parameters by 

n = 2Frf3h ls . (12) 

Thus larger values of k are accompanied by higher rates 
of decoherence through photon scattering, and this limits 
the time over which one can observe unitary evolution of 
the system. 

The kicking term in the Hamiltonian of Eq. (1) can 
be implemented by a train of magnetic field pulses sepa- 
rated by a time r. This results in a Zeeman interaction of 



the form g^sB ■ FJ2 n ^ — nT ) • The finite bandwidth 
limitations of magnetic coils and drivers prevent the ap- 
plication of true (5-kicks, but it is not difficult in practice 
to keep the kick duration T much shorter than the time r 
between kicks so that the 5-kick approximation remains 
valid. The angle of the turn per kick depends on the 
Larmor frequency of the applied magnetic field, 

P = n L r (13) 

Hence, by adjusting the laser intensity, detuning, Larmor 
frequency, kick spacing and duration, one can explore a 
whole range of kicked top parameters k and p. 

To observe quantum dynamics in different regions of 
the classical phase space, we start with initial spin coher- 
ent states |22 , which are rotations of the state \j, m = j) 
having maximum projection along the z-axis, 

\e,4>) = e mJ^<t>-Jyc°s4>]\j^ m = jy ( 14 ) 

The expectation value of the spin in this state is given 
by (J) = (jsmOcoscj), jsmBsmtp, jcosd). In the qubit de- 
scription, these states are separable with zero entangle- 
ment between qubits. In our atomic system, the total 
angular momentum F is the sum of a large nuclear spin 
(I = 7/2 for Cesium) and a valence electron spin-1/2 
system (qubit) whose reduced state can be described as 
in Eq. (7). Starting with initial spin coherent states 
and restricting ourselves to the manifold of maximum 
F (F — 4 for Cesium), the quantity S in Eq. (9) is 
then a measure of the entanglement between electron and 
nuclear spin, which can be experimentally measured by 
monitoring the mean collective angular momentum (F) 
through, e. g., Faraday spectroscopy [15] ■ Alternatively 
one can perform complete quantum state reconstruction 
of the overall electron+nuclear spin state [IS] . This opens 
up the possibility of computing other entanglement mea- 
sures such as negativity, and of monitoring decoherence 
by calculating the overall state purity. 

IV. DYNAMICAL ENTANGLEMENT AND 
CHAOS 

We analyze here the 1 : k entanglement dynamics of 
the kicked top described by the Hamiltonian in Eq. (1). 
We pick the magnitude of the angular momentum j to be 
4 to correspond to the F — 4 hyperfine manifold in the 
Cesium ground state. This puts the system far from the 
semi-classical regime of large j (j > 100) that has been 
studied in previous work on signatures of chaos [BJ \7\. 
Even in this deeply quantum regime with j = 4, we can 
clearly identify the effect of chaos in the entanglement 
dynamics. 

For k = 3,p = 7r/2, we pick an initial spin coher- 
ent state centered on a regular island of Fig. 1(b) with 
9 = 2.25,<f> = 2.5. The resulting evolution of the en- 
tanglement measure S shows quasi-periodic behaviour 
[Fig. 2(a)] with collapses and revivals in the entangle- 
ment. In contrast, for an initial state centered in the 
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FIG. 2: Evolution of the linear entropy S, (a) - (c) and the 
negativity N, (d) - (f) for initial spin coherent states \ 9, <j>) cen- 
tered in regular versus chaotic regimes of the classical phase 
space. See text for details. 



chaotic sea of Fig. 1(b) with 9 — 2.25,0 = 1.1, the 
quasiperiodic behaviour disappears [Fig. 2(b)] and the 
evolution is irregular. For this same initial state with 
9 = 2.25, <fi = 1.1, if we change n to 1, the classical 
dynamics becomes regular again [Fig. 1(a)], and corre- 
spondingly, quasiperiodic motion is recovered in the en- 
tanglement dynamics in Fig. 2(c). For both the regu- 
lar and chaotic dynamics, entanglement at first increases 
as the initially localized state starts spreading over the 
phase space. However, for the states initialized in regu- 
lar regimes, the dynamics causes the state to rclocalizc 
periodically, causing a reduction in entanglement. We 
explore this periodicity in more detail in the following 
section. 

We also compute a different measure of entanglement 
- namely the negativity defined as |23) 



N 



\Pt\ 



(15) 



where px is the partial transpose of the overall system 
state p, and the trace norm is defined to be 



\pt\\ = Tr 



PtPt 



(16) 



Entangled states have a non-zero negativity, but the con- 
verse is not necessarily true - i.e., states with zero neg- 
ativity may nevertheless be entangled. Whereas entropy 
is a measure of entanglement only for pure states, nega- 
tivity can be used as an entanglement measure for mixed 
states. We can thus use this measure later in Section V 
to understand the entanglement dynamics when we in- 
clude photon scattering in our model and the initial pure 
state becomes mixed due to decoherence. The evolution 
of the negativity for the three initial states considered 
in Fig. 2(a)-(c) is shown in Fig. 2(d)-(f). The qualita- 
tive behavior exactly matches the evolution of the linear 
entropy. We can identify collapses and revivals in the dy- 
namics for initial states in a regular regime and irregular 
evolution for states initially in the chaotic sea. 

The three initial states analyzed above are represen- 
tative of the general behavior of entanglement in regu- 
lar versus chaotic regions. The entanglement dynamics 
changes more or less smoothly from regular quasiperiodic 
evolution to irregular evolution as we scan through ini- 
tial spin coherent states from regular to chaotic regions. 
To illustrate this, Fig. 3(a) shows the average entangle- 
ment over 600 kicks as a function of initial conditions 
(9, <j>) scanned along the line 9 — 2.25. The quasiperiodic 
behavior for initial states in the regular islands leads to 
lower average entanglement as compared to initial states 
in the chaotic sea. In Fig. 3(a), the regular islands can 
be clearly identified by the dips in the average entangle- 
ment. Conversely, in order to confirm that the entangle- 
ment behaviour is connected to the level of chaoticity of 
the classical map, we scan through the chaoticity param- 
eter k for a fixed initial condition 9 = 2.25, cf> = 1.1. The 
time average entanglement as a function of n is shown 
in Fig. 3(b). As the chaoticity parameter is increased, 
the average entanglement increases. The dips in average 
entanglement reflect the occasional appearance of peri- 
odic orbits as fixed points in the classical phase space 
bifurcate in the approach to global chaos. 



V. ANALYSIS OF FLOQUET EIGENSTATES 

In previous work [B], we had identified similar signa- 
tures of chaos in the 2 : k entanglement dynamics of any 
two qubits with the remaining qubits for a collection of 
N > 50 qubits. Here, we have shown that these features 
are generic and can also be observed in the 1 : k entangle- 
ment between a single qubit and the rest. Furthermore, 
these signatures are surprisingly persistent even for the 
case j = 4 corresponding to only 8 qubits, which is far 
from the semi-classical regime. Ref. [7] explained these 
universal signatures of chaos by showing that in a semi- 
classical regime, the eigenstates \u n ), of the kick-to-kick 
Floquet evolution operator (Eq. 3) can be classified as 
'regular' or 'chaotic'. An initial state centered in a regu- 
lar island has support on a few almost degenerate regular 
eigenstates that give rise to quasiperiodic motion involv- 
ing a few regularly spaced cigcnfrcqucncies, whereas a 
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FIG. 3: Time averaged entanglement measured by both en- 
tropy (stars) and negativity (plusses) reflects regular and 
chaotic classical structures as (a) the location of the initial 
spin coherent state is varied along the line 6 = 2.25, < (j> < 
2tt for constant k = 3,p = n/2 and (b) k is varied keeping 
the initial state constant at \6 — 2.25, = 1.1). The average 
is taken over 600 kicks. For further details, see text. 

state initially localized in a chaotic sea can be decom- 
posed into a number of chaotic eigenstates with a broad 
spectrum of incommensurate frequencies contributing to 
the dynamics. 

Here, we test the above argument in the quantum 
regime. Figure 4 shows the eigenstate decomposition 
/ = | (u n \6, 4>) | 2 of the three initial states whose dynamics 
are shown in Fig. 2. Each eigenstate is labelled by the 
corresponding eigenphase ui n , 

U\u n )=e i ^\u n ), (17) 

where U is the Floquet evolution operator (Eq. 3). The 
regular initial state of Fig. 2 (a) and (d) has support on 
a few regular almost degenerate eigenstates [Fig. 4(a)], 
while the chaotic initial state of Fig. 2 (b) and (e) has 
support on eigenstates with a broader spectrum of ir- 
regularly space frequencies [Fig. 4(b)]. For the regular 
motion, the finite number of frequencies present in the 
power spectrum of S(t) can be identified as sums of dif- 
ferences between eigenfrequencies of the regular eigen- 
states on which the initial state has support [7]. In the 
chaotic case, the large number of distinct eigenfrequency- 



difference sums give rise to a broader power spectrum. 
The regular state in Fig. 4(c) also shows support mostly 
on a single eigenstate and some support on almost de- 
generate eigenstate pairs. Although the spectrum looks 
similar to the chaotic state of Fig. 4(b), the larger num- 
ber and irregular spread of eigenfrequencies in the chaotic 
state is sufficient to make the evolution much more irreg- 
ular than the regular case. Upon closer examination, we 
find that unlike the chaotic case, the eigenfrequency spac- 
ing of the dominant states in Fig. 4(c) is quite regular. 
This regular spacing and the larger support on a single 
eigenstate leads to regular entanglement dynamics, with 
periodic revivals or 'rephasing' of the evolutions of the 
different component eigenfrequencies, as seen in Fig. 2. 
Furthermore, support on the almost degenerate eigen- 
state pairs lead to long term periodic behaviour, which is 
missing in the chaotic case. The irregular oscillations on 
a fast timescale are due to the small but non-zero support 
on the remaining eigenstates. 

Due to the small value of j of 4, there are only a 
small total number of regular and chaotic eigenstates, 
so the chaotic power spectrum will not be completely 
flat and the dynamics will show some quasi-periodic be- 
haviour. Furthermore, due to the mixed nature of the 
phase space, some eigenstates have overlap on regular as 
well as chaotic regions and can partially contribute to 
initial states in both regular and chaotic regions. Thus 
the differences between regular and chaotic regimes are 
not as clearly delineated in this deeply quantum regime, 
compared to a semiclassical regime. Nevertheless, the 
regular and chaotic regions of the classical phase space 
can be identified by the support of the initial state on 
the Floquet eigenstates. To verify this, in Fig. 5 we plot 
the quantity 

s = ^|( U „|M>l (18) 

n 

as a function of initial conditions (9, <$) along the line 
9 = 2.25, for n = 3,p = 7r/2. This quantity roughly 
measures the number of eigenstates on which the initial 
state has support. The dips in s identify the regular 
regions in the mixed classical phase space and correlate 
perfectly with the time averaged entanglement. 

To further confirm the correspondence between the 
classical phase space and the quantum eigenstates, Fig- 
ure 6 shows the Husimi quasiprobability distribution of 
the eigenstats with largest overlap with each of the three 
initial states considered in Fig 2. The Husimi distri- 
butions show overlap of a state \u n ) with spin coherent 
states, 

P(9,cf,)= 2 ^\(9,cb\u n )\ 2 . (19) 

The distributions clearly show that the probability is con- 
centrated on the classical phase space structures shown 
in Fig 1. For regular eigenstates they are concentrated 
along periodic orbits whereas for chaotic eigenstates they 
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FIG. 4: Overlap / = \(u n \9, 4>)\ 2 of each of the 3 initial spin 
coherent states \0,<j>) in Fig. 2 with Floquet eigenstates \u n ) 
corresponding to eigenfrequencies ui n . Initial states in a reg- 
ular region (a), (c) have support on a few regularly spaced or 
degenerate eigenstates, while the state in the chaotic sea (b), 
has support on a larger number of chaotic eigenstates with 
irregularly spaced eigenfrequencies. 



are delocalized in the chaotic sea. Thus, although quan- 
tum chaos typically deals with semi-classical techniques 
to explore quantum-classical correspondence, our theo- 
retical analysis provides evidence that very clear signa- 
tures of classical chaos can be present in regimes far from 
the classical or even semi-classical regimes. 

A second signature of classical chaos predicted by the 
semiclassical analysis is that the initial increase in en- 
tanglement is faster for the state starting in the chaotic 
sea than for those starting in regular regions [7J. In the 
chaotic regime the initial rate of increase was predicted 
to be exponential due to the contribution of a large num- 
ber of incommensurate eigenfrequencies to the dynamics. 
In the quantum system considered here, due to the small 
size of the Hilbert space, there are not enough eigenfre- 
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FIG. 5: Support on the Floquet eigenstates, |lt n ), s = 
5~Jn |(u n |0, 4>)\. The initial spin coherent state \0, 4>) is varied 
along the line 6 = 2.25, < <f> < 2tt with fixed k = 3, p = n/2. 
s is normalized with respect to its value of so = 3, corre- 
sponding to the case of equal support over all eigenstates. 
The dips correspond to the location of the regular islands in 
the classical phase space of Fig. 1. 



quencies to give rise to an initial exponential increase 
in entanglement. Furthermore, the initial spin coherent 
states are not well localized in the phase space and there- 
fore tend to have support on both regular and chaotic 
regions. This in turn tends to wash out differences in 
the initial rate of entanglement generation. In fact, we 
find no significant difference in the initial increase in en- 
tanglement for states in regular versus chaotic regions. 
Our studies thus not only serve to identify the corre- 
spondence between quantum and classical descriptions, 
but also highlight the differences inherent in quantum 
and classical systems. 



VI. DYNAMICAL ENTANGLEMENT AND 
DECOHERENCE 

In order to accurately model any physical system such 
as, e. g., Cesium atomic spins, we must include the effects 
of decoherence due to interaction with the environment. 
We perform simulations for the specific case of the cold 
atom system including photon scattering which results 
in non-unitary evolution. We model the evolution of the 
density operator p of the system using a master equation 
of the form 



dp 
At 



UH, P ] + Y t WpD\-\{D\D q ,p} 



(20) 



e q ■ d ge describe spontaneous 



The jump operators D q 

emission from the excited states [53]. The excited states 
are adiabatically eliminated by assuming that the ground 
state populations change slowly compared to the excited 
states and the coherences. The excited state populations 
are thus effectively slaved to the ground state populations 
to give a net evolution only in terms of the ground state 
manifold F = 4. Furthermore, any population which 
is pumped into the F = 3 ground state is treated as 
loss and ignored in the subsequent dynamics. This is 
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j=4, k=1 ,p=jt/2,8=2.25,cp=1 .1 




FIG. 6: Husimi distribution of the Floquet eigenstate having 
the largest overlap with each of the three initial states (a)-(c) 
considered in Fig. 2. Solid lines are contours of equiproba- 
bility distribution, and darker regions indicate larger proba- 
bility. Comparison with Fig. 1 clearly shows concentration of 
the distributions along the classical phase space structures. 



possible if atoms in the F = 3 manifold are not re-excited 
and therefore do not participate further in the observed 
dynamics. This model has provided accurate simulations 
of the physical system in recent experiments [TBI 1X81 IH)] . 



A. Effect of decoherence on entanglement 

Figure 7 shows the simulated evolution (Eq. 20) of the 
entropy S(t) including decoherence for the three initial 
states analyzed in Fig. 2. Comparing the unitary evolu- 
tion of Fig. 2(a)-(c) to the non-unitary evolution due to 
decoherence in Fig. 7 (a)-(c), we see that the qualitative 
features are the same. The initial states in the regular 
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FIG. 7: Evolution of the linear entropy S, (a) - (c) and the 
negativity N, (d) - (f) in the presence of decoherence (photon 
scattering) for initial spin coherent states \6,(j>) centered in 
regular versus chaotic regimes of the classical phase space. 
For comparison, the 3 initial states considered are picked to 
be the same as those in Fig. 2. See text for details. 



regimes give rise to quasiperiodic motion, which is miss- 
ing in the chaotic regime. Thus the signatures of chaos 
discussed in the previous section are also evident in the 
open system dynamics. However, the total time of evo- 
lution is limited by photon scattering, which causes the 
system to get entangled with the environment. The lin- 
ear entropy no longer quantifies the 1 : k entanglement of 
a single qubit (electron) with the rest (nucleus). Instead 
it quantifies the total entanglement with both the system 
and the environment. As the system becomes more en- 
tangled with the environment, there is an overall increase 
in entropy towards its maximum value of 0.5. 

In order to better understand the effect of decoherence 
on the signatures of chaos in the entanglement dynam- 
ics, we compute the negativity as given by Eq. (15). Un- 
like the entropy S, the negativity N does give us a valid 
measure of 1 : k entanglement for an overall mixed state 
resulting from tracing out the environment. Fig. 7(d)- 
(f) shows the evolution of the negativity for the initial 
states considered in Fig. 2. We see that there is an 
overall decrease in the negativity but the quasiperiodic 
behavior in the regular regimes and the irregular dynam- 
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ics in the chaotic regime exactly match the evolution of 
the entropy S in Fig. 7(a)-(c). This shows that the to- 
tal spin state collectively gets more entangled with the 
environment, and accordingly the total 1 : k entangle- 
ment within the system decreases. However, because de- 
coherence acts collectively on the total spin system, the 
signatures of chaos caused by the interactions within the 
system are not completely destroyed and persist well be- 
yond the scattering time. 

Although negativity has the advantage that it is a valid 
measure of entanglement for mixed states, it presents a 
disadvantage for experiments when compared to the lin- 
ear entropy. This is because, in order to study negativity 
dynamics in an experiment, the total (nuclear + elec- 
tron) state must be tomographically reconstructed. This 
is experimentally more challenging than measuring the 
linear entropy via measurements of (F a ) (Eq. 9). Hence, 
in actual experiments, it may be more convenient only to 
measure the linear entropy S. Our comparison of the neg- 
ativity and the linear entropy (Fig. 7) then clarifies that 
the decoherence only causes an overall decrease in the 
entanglement between electron and nuclear spin. Thus 
the structures experimentally observable in the linear en- 
tropy dynamics would correspond to valid signatures of 
chaos in the entanglement between electron and nuclear 
spin. Furthermore the overall effect of the environment 
averages out in the time averaged entanglement dynam- 
ics, leading to curves very similar to the unitary evolution 
in Fig. 3. 



B. Effect of chaos on decoherence rate 

In addition to examining the entanglement within the 
spin system, we can also analyze the entanglement of 
the total spin state with the environment due to photon 
scattering. We set the intensity, detuning and Larmor 
frequency to obtain k = 3 and p = ir/2. We then compare 
the decoherence rates of a spin coherent state centered on 
the regular island at 9 = 2.25, cj) = 2.5 and a spin coherent 
state centered in the chaotic sea at = 2.25 and cf) = 1.1. 
Since the experimental parameters are the same for both 
cases, the scattering rate 7 S is also the same. However, 
despite the scattering rate being constant, we find that 
the decoherence rate is faster for the state initiated in the 
chaotic sea relative to the state starting on the regular 
island. In Fig. 8, we plot the purity of the overall state 
as a function of time. The purity decay is a measure of 
the decoherence and shows that the decoherence rate is 
faster for a state in the chaotic sea than for a state on a 
regular island. We have observed this increased rate of 
decoherence for states in a chaotic sea for other values of 
k and p as well, indicating that chaos generally acts to 
increase the decoherence rate of the system. 

The increase in the decoherence rate in the chaotic 
regime can be understood by examining the Lindblad 
terms in the master equation of Eq. (20) more carefully. 
The populations in the excited states which determine 
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FIG. 8: Purity decay (decoherence) of the overall spin state 
for different initial conditions with fixed k, = 3,p = 7r/2. The 
decoherence rate is faster for an initial spin coherent state 
centered on a regular island with \6 = 2.25, cj> = 2.5} (dashed) 
compared to an initial spin coherent state centered in the 
chaotic sea with \8 = 2.25,0 = 1.1) (solid), although the 
photon scattering rate is the same for both cases. 



the spontaneous emission rate, depend on the ground 
state populations. Figure 9 shows the Husimi distribu- 
tions after 50 kicks for the two initial spin coherent states 
centered in the regular island and chaotic sea. Whereas, 
the initial state on a regular island remains close to a 
more robust spin coherent state, for the initial state in 
the chaotic sea, the dynamics causes the population dis- 
tribution to quickly lose any symmetry and diverge from 
a coherent state. The spontaneous emission for such a 
state leads to a net larger decoherence rate. 

Our results support the arguments originally formu- 
lated by Zurek and Paz [HI [12] claiming that chaos can 
enhance the rate of decoherence. The relevance of our 
work is that we consider a realistic and accurate decoher- 
ence model for an actual physical system very deeply in 
the quantum regime and verify the predictions obtained 
from semiclassical theories and general decoherence mod- 
els. 



VII. CONCLUSION 

In summary, we have presented an analysis of entan- 
glement and decoherence in the quantum kicked top in a 
regime of relevance to a proposed experiment with cold 
Cesium atoms. 1 : k entanglement can be easily cal- 
culated from the mean angular momentum vector, and 
signatures of chaos are evident in the dynamics of en- 
tanglement in a quantum regime. The dynamics can be 
understood by examining the decomposition of the initial 
state into regular and chaotic cigenstates of the Floquct 
operator. 

Our simulations of a realistic system of cold atoms 
undergoing decoherence due to photon scattering reveal 
that decoherence reduces the 1 : k entanglement but does 
not erase the signatures of chaos. Furthermore, we have 
shown that chaos can enhance the overall decoherence 
rate, or entanglement of the system with the environ- 
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ment. This shows that whereas chaos can be helpful for 
generating entanglement within the system, it can have 
negative effects by rapidly causing decoherence. 



We have thus shown a means for performing the first 
experimental studies of entanglement in a chaotic sys- 
tem, identified signatures of chaos in entanglement in a 
quantum regime that previous theoretical studies have 
not investigated, and extended the previous theoretical 
analysis to this quantum regime. Furthermore, we have 
demonstrated the collective effect of decoherence on the 
entanglement, and identified and explained signatures of 
chaos in the decoherence rate using a realistic model of 
decoherence. An effort to implement the experiments de- 
cribed here is currently underway. 



FIG. 9: Husimi distributions after 50 kicks for the initial 
states considered in Fig. 8, with k = 3,p = n/2. Solid 
lines are contours of equiprobability distribution, and darker 
regions indicate larger probability. The initial spin coherent 
state centered on a regular island (a) with \0 = 2.25.0 = 2.5) 
remains closer to a robust spin coherent state, whereas the 
initial spin coherent state centered in the chaotic sea with 
\9 — 2.25.0 = 1.1} spreads quickly, leading to a delocalized 
state with net larger spontaneous emission and hence larger 
decoherence. 
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